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Abstract 



It is shown, that the energy density of coherent axion field oscillations in the cos- 
■ mology of standard invisible axion should be distributed in the Universe in the form 

of archioles, being nonlinear inhomogeneous structure, reflecting the large scale distri- 
ct bution of Brownian structure of axion strings in the very early Universe. Spectrum of 
inhomogeneities, generated by archioles, is obtained and their effects in the spectrum and 
quadrupole anisotropy of relic radiation are considered. The axionic-string-decay-mo del- 
independent restriction on the scale of axion interaction is obtained. 



1 Introduction 



The modern theory of large-scale structure of the Universe is based on the assumption that this 
structure is formed as the result of development of gravitational instability from small initial 
perturbations of density or gravitational potential. As a rule, these perturbations are gaussian, 
but some version of nongaussian perturbations have also been considered. 

In this paper we analyze the problem, inherent to practically all the cosmological cold dark 
matter models of invisible axion, that concerns primordial inhomogeneity in the distribution 
of the energy of coherent oscillations of the axion field. This problem, referred to as the 
problem of archioles, invokes nongaussian component in the initial perturbations for axionic 
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cold dark matter. Archioles are the formation that represents a replica of the percolation 
Brownian vacuum structure of axionic walls bounded by strings, which is fixed in the strongly 
inhomogeneous primeval distribution of cold dark matter. They can give rise to interesting 
alternative scenarios of structure formation that relate the mechanism responsible for structure 
formation to inhomogeneities of the type of topological defects. 

The analysis of observable effects associated with archioles leads to a new model-independent 
constraint on the mass of invisible axion. Such analysis is also very useful for further develop- 
ment of full cosmological theories, based on the model of horizontal unification (MHU), which 
has been proposed in Ref. |IJ as the minimal phenomenology of everything, including the physics 
of inflation, bariosynthesis, and dark matter. In particular, the combination of archioles effect 
with the consideration of nonthermal symmetry restoration in the horizontal phase transitions 
on inflation stage, puts the upper limit on the scale of family symmetry breaking (the main 
parameter of MHU) and consequently severely reduces the set of possible realizations of dark 
matter scenarios in this model. 



2 Formation of the archiole structure in the early Universe 

In the standard invisible axion scenario Ref. [[| the breaking of the Peccei- Quinn symmetry is 
induced by the complex SU(3) (g) SU(2) ® U{\) - singlet Higgs field with a "Mexican hat" 
potential 

V{<P) = \^ + ct>-Fl) 2 (1) 

Such field can be represented as 4> = F a exp(i9), where 9 = a/F a and a is the angular Gold- 
stone mode-axion. QCD instanton effects remove the vacuum degeneracy and induce effective 
potential for 9 

V{9) = k\(l - cos(9N)) (2) 

Below we will simply assume for standard axion that N — 1. In the context of standard big 
bang scenario it is usually assumed that the phase transition with U(l) - symmetry breaking 
occurs when the Universe cools below the temperature T = F a . Thus, in the standard case the 
crucial assumption is that from the moment of the PQ phase transition and all the way down to 
the temperatures T = Aqcd, the bottom of the potential Eq. (P is exactly flat and there is no 
preferred value of a during this period (the term given by Eq. (|2|) vanishes). Consequently, at the 
moment of the QCD phase transition, when the instanton effects remove vacuum degeneracy, a 
rolls to the minimum and starts coherent oscillations (CO) about it with energy density Ref. [FJ] 

p „ M = 19 . 57 (^)(£) 3 r^ (3) 

The coherent axion field oscillations turn on at the moment t « 8.8 ■ 10~ 7 s. 

It is generally assumed, that PQ transition takes place after inflation and the axion field 
starts oscillations with different phase in each region causally connected at T = F a , so one 
has the average over all the values to obtain the modern axion density. Thus in the standard 
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cosmology of invisible axion, it is usually assumed that the energy density of coherent oscil- 
lations is distributed uniformly and that it corresponds to the averaged phase value of 9 = 1 
(p a = p{9)). However, the local value of the energy density of coherent oscillations depends on 
the local phase 9 that determines the local amplitude of these coherent oscillations. It was first 
found in Ref. 0, that the initial large-scale (LS) inhomogeneity of the distribution of 9 must 
be reflected in the distribution of the energy density of coherent oscillations of the axion field. 
Such LS modulation of the distribution of the phase 9 and consequently of the energy density 
of CO appears when we take into account the vacuum structures leading to the system of axion 
topological defects. 

As soon as the temperature of Universe becomes less then F a , the field <ft acquires the 
vacuum expectation value (VEV) (0) = -F a exp (i9), where 9 varies smoothly at the scale F~ l . 
The existence of noncontractable closed loops that change the phase by 27m leads to emergence 
of axion strings. These strings can be infinite or closed. The numerical simulation of global 
string formation Ref. M revealed that about 80% of the length of strings corresponds to infinite 
Brownian lines. The remaining 20% of this length is contributed by closed loops. Infinite 
strings form a random Brownian network with the step L(t) ~ t. After string formation when 
the temperature becomes as low as T « A-qcd, the Eq. ([|) makes a significant contribution 
to the total potential so that the minimum of energy corresponds to a vacuum with 9 = 2nk, 
where k is an integer - for example, k = 0. However, the vacuum value of the 9 cannot 
be zero everywhere, since the phase must change by A9 = 2tt upon a loop around a string. 
Hence, we come from the vacuum with 9 = to the vacuum with 9 = 2tt as the result of such 
circumvention. The vacuum value of 9 is fixed at all points with the exception of the point 
9 = 7i. At this point, a transition from one vacuum to another occurs, and the vacuum axion 
wall is formed simultaneously with CO turning on. The width of such wall bounded by string 
is S = m^ 1 . Thus, the initial value of 9 must be close to tt near the wall, and the amplitude 
of CO in Eq. is determined by the difference of the initial local phase 9(x) and the vacuum 
value, which is different from the one of the true vacuum only in a narrow region within the 
wall of thickness 5 = m~ l . Therefore in this region we can write 9(x) = tt — e(x), where Ref. || 
e(x) = 2 tan _1 (exp(m a :r)) and x = m" 1 . Thereby the energy density of CO in such regions is 
given by 

P A ~ vr 2 pa (4) 

And so we obtain, that the distribution of CO of axion field is modulated by nonlinear inho- 
mogeneities in which relative density contrasts are 5p/p > 1. Such inhomogeneities were called 
archioles. In the other words archioles are a formation that represents a replica of the perco- 
lational Brownian vacuum structure of axionic walls bounded by strings an which is fixed in 
the strongly inhomogeneous initial distribution of axionic CDM. The scale of this modulation 
of density distribution exceeds the cosmological horizon because of the presence of 80% infinite 
component in the structure of axionic wall bounded by strings system. The superweakness of 
the axion field selfinteraction results in the separation of archioles and of the vacuum structure 
of axionic walls - bounded - by - strings. So these two structures evolve independently. The 
structure of walls bounded by strings disappears rapidly due to disintegration into separate 
fragments and further axion emission. The structure of archioles remains frozen at the RD 
stage. On the large scales, the structure of archioles is an initially nonlinear formation, a 
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Brownian network of quasi - one - dimensional filaments of dustlike matter with the step 



L A (t) = Xt (5) 

(where A = 1). At the moment of creation t, the linear density of this quasilinear filamentary 
formations given by 

Ha = 7!~ 2 pj5 (6) 

In accordance with this, the cosmological evolution of archioles in the expanding Universe is 
reduced to the extension of lines along only one direction. 

We have studied in Ref. 0] the spectrum of inhomogeneities that the density develops in 
response to the large-scale Brownian modulation of the distribution of CO of axion field. Density 
perturbations, associated with Brownian network of archioles, may be described in the terms 
of a two - point autocorrelation function Ref. |S|. To obtain such autocorrelation function, 
it is necessary to perform averaging of energy density of infinite Brownian lines over all lines 
and over the Winner measure, which corresponds to the position along of Brownian line (see 
Ref.0). 

The two - point autocorrelation function in the Fourier representation has the form 

(^-(k) 5 -^-{k')) = l2pAfi A k' 2 5(k + k')r l r 2 t A G 2 (7) 
Po Po 

where p is background density, Jmd = 3/(327r) for dustlike stage, fnc = (67r) _1 for RD stage, 
G is the gravitational constant, pa is the total energy density of the Brownian lines. The mean- 
square fluctuation of the mass is given by 

= UpAPAt-'r'G'kt 4 (8) 

3 Cosmological impact of archioles 

Let us consider a region characterized at instant t by a size I and a density fluctuation A. For 
anisotropy of relic radiation we then obtain 




ST 
~T 




(9) 



If Z = t, we have \ST/T\ = |A|; that is, the anisotropy of relic radiation is equal to the density 
contrast calculated at the instant when the size of the region is equal to the size of the horizon 
(Sachs - Wolf effect). To estimate the quadrupole anisotropy that is induced in relic radiation 
by the structure of archioles, we must find the amplitude of perturbations on the scale of the 
modern horizon 
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Thus Sachs- Wolf quadrupole anisotropy of relic radiation induced by archioles will be 

According to COBE data (see for examle Ref. @), the measured quadrupole anisotropy of relic 
radiation is at the level of 

Y « 5 • 1(r6 (12) 

If we take into account the uncertainties of our consideration such as the uncertainties in 
correlation length scale of Brownian network (A ~ 1 4- 13) and in temperature dependence of 
axion mass, we can obtain a constraint on the scale of symmetry breaking in the model of an 
invisible axion 

F a < 1.5- 10 10 GeF 4 4 ■ 10 9 GeV; m a > AlOfxeV 4 1500/iel/ (13) 

This upper limit for F a is close to the strongest upper limits in Refs. JF| [8|, obtained by 
comparing the density of axions from decays of axionic strings with the critical density, but has 
an essentially different character. 

The point is that the density of axions formed in decays of axionic strings depends critically 
on the assumption about the spectrum of such axions (see Refs. ]7|, [SJ) and on the model of 
axion radiation from the strings (see Ref. ||). For example, Davis Ref. @ assumed that radiated 
axions have a maximum wavelength of uj(t) = t^ 1 while Harari and Sikivie Ref. M have argued 
that the motion of global strings was overdamped, leading to an axion spectrum emitted from 
infinite strings or loops with a flat frequency spectrum oc k . This leads to an uncertainty 
factor of ~ 100 in the estimate of the density of axions from strings and to the corresponding 
uncertainty in the estimated upper limit on F a 

F a <2 - 10 1C VGW; m a > 300/^eV (14) 

Here, ? = 1 for the spectrum from Davis Ref. |7| , and <j ~ 70 for the spectrum from Harari and 
Sikivie Ref. §. 

In their treatment of axion radiation from global strings, Battye and Shellard Ref. found 
that the dominant source of axion radiation are string loops rather than long strings, contrary 
to what was assumed by Davis Ref. 0. This leads to the estimations 

F a < 6 • !0 10 GeV 4 1.9 • !0 n GeV; m a > 3lfxeV 4 lOOfieV (15) 

Arguments that lead to the constraint Eq. ( |T3"D are free from these uncertainties, since they 
have a global string decay model - independent character. 

At the smallest scales, corresponding to the horizon in the period t, evolution of archioles 
just in the beginning of axionic CDM dominancy in the Universe (at redshifts zmd — 4 • 10 4 ) 
should lead to formation of the smallest gravitationally bound axionic objects with the minimal 
mass M ~ p a t 3 ~ 10~ 6 M Q and of typical minimal size t(l + za) /(l + %md) — 10 13 cm. One can 
expect the mass distribution of axionic objects at small scale to peak around the minimal mass, 
so that the existence of halo objects with the mass (1O _6 M ~ 1O _1 M ) and size 10 13 4- 10 15 cm 
is rather probable, what may have interesting application to the theoretical interpretation of 
MACHOs microlensing events. 
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4 Physical impact of archioles 



The inclusion of obtained restriction into the full cosmoparticle analysis provides detailed quan- 
titative definition of the cosmological scenario, based on the respective particle physics model. 
Consider, for example, a simple variant of gauge theory of broken family symmetry (TBFS) 
Ref. [[UJ, which is based on the standard model of electroweak interactions and QCD, sup- 
plemented by spontaneously broken local SU(3)h symmetry for quark-lepton families. This 
theory provides natural inclusion of Peccei-Quinn symmetry U(1)h = U(1)pq, being associ- 
ated with heavy "horizontal" Higgs fields and it gives natural solution for QCD CP - violation 
problem. The global U(1)h symmetry breaking results in the existence of axion-like Goldstone 
boson - a. TBFS turns to be a simplest version of the unified theoretical physical quantitative 
description of all main types of dark matter (HDM-massive neutrinos, axionic CDM and UDM 
in the form of unstable neutrinos Refs. [|H], |Tl|) and the dominant form of the dark matter is 
basically determined by the scale of the "horizontal" symmetry breaking Vh, being the new 
fundamental energy scale of the particle theory. For given value of Vh the model defines the 
relative contribution of hot, cold and unstable dark matter into the total density. Since in 
the TBFS the scale of horizontal symmetry breaking Vh is associated with F a , we have, from 
Eq. (0), the same upper limit on Vh- However, this limit assumes, that the considered infla- 
tionary model permits topological defects and hence archioles formation due to the sufficiently 
high reheating temperature Trh > Vh- In the inflationary model, which occurs in TBFS, we 
can achieve Trh ~ 10 10 GeV. The "horizontal" phase transitions on inflationary stage lead 
to the appearance of a characteristic spikes in the spectrum of initial density perturbations. 
These spike-like perturbations, on scales that cross the horizon (60-1) e- folds before the end of 
inflation reenter the horizon during the radiation or dust like era and could in principle collapse 
to form primordial black holes. The minimal interaction of "horizontal" scalars of TBFS 
£W, £( 2 ) with inflaton allows us to include them in the effective inflationary potential Ref. [Tj|: 

e<°>, e«, e (2) ) = -^0 2 + V - E ^ (e«) 2 + e4" (^) 4 + E v h 2 (^f ( 16 ) 

1 4 i=0 1 i=0 4 i=0 1 

The analysis of processes of primordial black holes formation from density fluctuations, which 
can be generated by "horizontal" phase transitions at the inflationary stage gives rise to an 
upper limit on the scale of horizontal symmetry breaking Ref. |12| . 

V H < 1.4 ■ 10 13 GeV (17) 

Therefore the range between the two upper limit Eq. (|IB|) and Eq. ( [T7| ) turns to be not 
closed, and the following values seem to be possible 

10 u GeV <V H < 10 13 GW (18) 

The indicated range corresponds to the case when all the horizontal phase transitions take 
place on the dust like stage and C2 Cmp;. In this case the inflationary field <fi oscillates with 



initial amplitude ~ mpi. According to Refs. [[L3, [L2| it means that any time the amplitude of 



the field becomes smaller then C2 <C mpi, the last (axion £^ 2 ^) phase transition with symmetry 
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breaking occurs, and topological defects are produced. Then the amplitude of the oscillating 
field becomes greater than C2 , and the symmetry is restored again. However, this regime 
does not continue too long. Within a few oscillations, quantum fluctuations of the field ^ will 

be generated with the dispersion ((^ 2 - ) ) ) — v^ 1 X^ 2 ln~ 2 1/V| Ref. [fl3f1 . For 



ml^u^X^X^ln-H/ul (19) 

these fluctuations will keep the symmetry restored. The symmetry breaking will be finally 

completed when (k^) ) will become small enough. Thus such phase transition leads to 
formation of topological defects and archioles without any need for high - temperature effects. 
Substituting the typical values for potential Eq. ( |TE| ) such as m| w 10~ 3 V^, X^ ~ 10~ 3 , — 
1CT 10 (see Ref. ||12|| ) we will obtain that the condition Eq. (|i~9"D means that for the scales 



V H < KT 3 m Pi (20) 

the phenomenon of non - thermal symmetry restoration takes place in simplest inflationary 
scenario based on TBFS. Owing to this phenomenon oscillations of the field ^ do not suppress 
the topological defects and archioles production for the range Eq. (pT8[). So the range Eq. (|1"8| ) 
turns to be closed by comparison of BBBR quadrupole anisotropy, induced by archioles, with 
the COBE data. As a result, the upper limit on the scale of horizontal symmetry breaking will 
be given by Eq. (|13|). 

Note, in conclusion, that the axion emission can influence the time scale and energetics of 
neutrino flux from collapsing stars. Analysis of this effect for SN1987A excludes the interval 
3-10 6 Ge^ < Vh < 3T0 9 CreV (see Ref. | JL4|| ) and establishes the lower limit on the high energetic 
branch of TBFS. Thus putting together all these limits we can extract narrow window in the 
high energetic branch of the so called model of horizontal unification (MHU) Ref. 

3 • 10 9 GeV <V H < 1.5 • 10 w GeV (21) 

On the base of this choice for the main parameter of MHU we can build a quantitatively 
definite dark matter scenario, which associates the formation of the large-scale structure in 
the Universe with a mixture of axions and massive neutrinos, since in this interval the total 
density equal to the critical one makes in the most cases the contribution of massive neutrinos 
necessary. 
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